On Spherically Symmetric Motions of the 
Atmosphere Surrounding a Planet Governed 
by the Compressible Euler Equations 



We consider spherically symmetric motions of inviscid compress- 
ible gas surrounding a solid ball under the gravity of the core. Equi- 
libria touch the vacuum with finite radii, and the linearized equation 
around one of the equilibria has time-periodic solutions. To justify the 
linearization, we should construct true solutions for which this time- 
periodic solution plus the equilibrium is the first approximation. But 
this leads us to difficulty caused by singularities at the free bound- 
ary touching the vacuum. We solve this problem by the Nash-Moser 
theorem. 
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1 Introduction 

We consider spherically symmetric motions of atmosphere governed by the 
compressible Euler equations: 
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Abstract 
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1 



and the boundary value condition 

pu\ r =R = 0. (2) 

Here p is the density, u the velocity, P the pressure. Rq{> 0) is the radius 
of the central solid ball, and g = G M , G being the gravitational con- 
stant, M the mass of the central ball. The self-gravity of the atmosphere is 
neglected. 

In this study we always assume that 

P = Ap\ (3) 

where A and 7 are positive constants, and we assume that 1 < 7 < 2. 

Equilibria of the problem are given by 

«r) = JM;-5)* <*^<*> 

(0 (R<r), 
where R is an arbitrary number such that R > Ro and 

Remark: The total mass M of the equilibrium is given by 



JR ^ R, 

M is an increasing function of R. Of course M — > as R — > R . But as 
R — > +00, we see 

*, a a r 22^1 J f+°° if 7 > 4/3 

M ->■ A-kAi r 7-1 dr = < ' 

Jr \M*(<oo) if 7 < 4/3, 

where 

= ij p 7 _! 

4-37 

Hence if 7 > 4/3 there is an equilibrium for any given total mass, but if 
7 < 4/3 the possible mass have the upper boud M*. 
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Let us fix one of these equilibria. We are interested in motions around this 
equilibrium. The major difficulty of the analysis comes from the free bound- 
ary touching the vacuum, which can move along time. So it is convenient to 
introduce the Lagrangian mass coordinate 



= 47r / p(t,r')i 



m = 47r / p(t, r )r dr , 
'Ro 

to fix the interval of independent variable to consider. Taking m as the 
independent variable instead of r, the equations turn out to be 

du A 2 dP g , n „,s 



where 



We note that 



/ 3 [ m dm\Vz 



dr dr 1 



dt dm Ai\r 2 p 

Let us take f = f(m) as the independent variable instead of m, where 
f{m) is the inverse function of 

m(r) = An I p{r')r' 2 dr' . 
JRo 

Then, since 

d Id / 2 dr\~ x _/r 2 dr\- 1 

dm Ai[f 2 pdf' ^ V dm) ^\f 2 df) ' 

we have a single second-order equation 

d 2 r lr^^/-/r^dr\--i\ + 9o = Q 
dt 2 pf 2 df\ \f 2 dfJ J r 2 

The variable f runs on the interval [Ro, R] and the boundary condition is 

r\f=R = Ro- 
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Without loss of generality, we can and shall asume that 

R = 1, g = — *— -, A — -, A ± = 1. 
7-1 7 



Keeping in mind that the equilibrium satisfies 

p dr r 2 

we have 



<9 2 r 1 r 2 <9 fp/^ /r 2 <9r\ -t\ \ 1 /l r 2v 



<9t 2 pr 2 dr\ V \r 2 dr / 7/ 7 — 1 Vr 2 r 4 , 
Introducing the unknown variable y for perturbation by 

r = r(l + y), (4) 

we can write the equation as 
where 

G(y,v) := l-(l + y)- 2 T(l+y + t;)^ = 7(3y + t;) + [i/,t;] 2 , 
lf(y) := (l + y) 2 -^-^ = 4y + [y] 2 

and we have used the abbreviations r, p, P for f, p, P. 

Notational Remark: Here and hereafter [X\ q denotes a convergent 
power series, or an analytic function given by the series, of the form Ylj> q a jX\ 
and [X, Y] q stands for a convergent double power series of the form J2j+k>q a jk^Y k . 

We are going to study the equation (5) on 1 < r < R with the boundary 
condition 

y\ r =i = 0. 
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Of course y and r— will be confined to 
or 







\y\ + 


dr 



< 1. 



2 Analysis of the linearized problem 
2.1 

The linearization of the equation (5) is clearly 

d 2 y 



dt 2 



+ Cy = 0, 



(6) 



where 



Cy ■= - 



1 d 

pr dr 



(1 l\d 2 y / 4/1 1\ 7 1 \ rfy 37 - 4 y 
r i?/ rfr 2 V r Vr i?/ 7 — 1 r 2 / dr 7 — 1 r 3 



Let us consider the eigenvalue problem 

Cy = \y. 

We introduce the independent variable z by 



(7) 



R 



z = 



R 



and the parameter N by 
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7-1 



iV 2 
— or 7 = 1 H . 

2 [ iV-2 



Then we can write 
R 3 Cy = - 



z d' 2 y f - Az dy 8 - N I 



1 — z dz 2 (\ — z) 2 dz 



1-zf 



(8) 



(9) 



(10) 



The variable z runs over the inetrval [0,1 — -|], the boundary z = 
corresponds to the free boundary touching the vacuum, and the boundary 
condition at z = 1 — is the Dirichlet condition y — 0. 



2.2 

Although the boundary z = 1 — ^ is regular, the boundary 2 = is singular. 
In order to analyze the singularity, we transform the equation Cy — Xy , 
which can be written as 



n, ± 1_1 

(1- :) 

to an equation of the formally self-adjoint form 



^ (N 1 4^ xdy 8-JV y 3 



This can be done by putting 



8-jV 
I 2 



— / 

p = z 2 (1 — z 

8 — iV N-2 4-Af 

9 = Z 2 (1 - Z) 2 . 



w — ^ . 
/i = Z 2 (1 — Z) 



Using the Liouville transformation, we convert the equation 



to the standard form 



de 

This can be done by putting 



d2,] + Qn = xr 3 ,, + f. 



v = M 1/4 y = Z 4 (1 - Z) 4 y, 

/ = p 1 /v 3 /V = ^(i-^)^/, 

and 



/i A\p /x / loVp /i / Ap\p \i J J 

- 1 ^ (^-l)(^-3) 7-2iV x 
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Putting 



£ 2 

we see the variable £ runs over the inetrval [0,£r]. Since z ~ — as £ — > 0, 

we see 

(iV-l)(iV-3) 1 

W ~ 4 e 

, (JV-l)(iV-3) 3 

as t — > 0. But 7 < 2 implies N > 4 and > -. Hence the 

1 4 4 

boundary £ = is of the limit point type. See, e.g., [4J, p. 159, Theorem X.10. 
The exceptional case 7 = 2 or N = 4 will be considered separately. Thus we 
have 

Proposition 1 The operator T ,V(T ) = C^°(0, T r] = —rj^ + Qr], in 
L 2 (0,^r) has the Friedrichs extension T, a self-adjoint operator whose spec- 
trum consists of simple eigenvalues X1R 3 < X2R 3 < ■ ■ ■ < A n i? 3 <•••—» +00. 

In other words, the operator Aq,V(Ao) = C£°(0, 1 — ^),A y = Cy in 

X := L 2 ((0, 1 — j^), fidz(= z~2^(l — z) 12 ^ dz)) , has the Friedrichs extension 
A, a self-adjoint operator with the eigenvalues (A n ) n . 
We note that the domain of A is 

V(A) = { y eX\3<p n eC™(0,l-±) 

such that (p n — > y in X and 0[</> m — <f> n ] — > as m, n — > 00, 
and Cy G X in distribution sense.} 

Here 



Of 



r 

^0 


d(f) 
dz 


2 A 

Z 2 


(1- 


-r 


1 

R 


z 




Jo 


1 


— z 


dz 



2 

ji(z)dz. 



Moreover we have 

Proposition 2IfN<8(or , ~f> 4/3 J ; the least eigenvalue Ai zs positive. 
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Proof: Suppose N < 8. Clearly y = 1 satisfies 

d dy 8 — N n-2 a-n 

-Tz p Tz + qy = q = — z 2 {1 ~ z) 2 ' 

Therefore the corresponding 771 (£) given by 

JV-l 

771=2; 4 (1 - Z) 

satisfies 

d 2 Vi „ . 8-N 1 



9-N 
4 



+ Qv — Q — z 4 (1 — 4 • 



de * 2 

It is easy to see 771 = -j^ = at £ = 0. Let 0i(£) be the eigenfunction of 
—d 2 /dC, 2 + Q associated with the least eigenvalue Ai. We can assume that 
<M0 > for < £ < U,MU) = 0, and < at £ = £ R . Then 
integrations by parts gives 



> 



> 0. 



Remark: When N — 8, 771 satisfies — 37^ + Q77 = 0, but does not satisfy 



(i 2 77 

the boundary condition 77^=^ = 0. Hence it is not an eigenfunction of zero 
eigenvalue, and Ai > even if N = 8. 
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2.3 



For the sake of convenience of the further analysis, let us rewrite the linear 
part C by introducing a new variable 

Clearly x = z+ [z] 2 and the change of variables z t-t xis analytic on < z < 1 
and its inverse x >->■ z is analytic on < x < n 2 /16. Since 

d I y/l ^ 

dz V z dx'' 

d 2 1 — z ( d 2 If [x~ 1 



1 + 



cb 2 z V cfe 2 2\ V z {\ — z)\fY^z ) dx, 

we can write 

where £i(5t) and ^o(^) are analytic on < x < n 2 /16. Putting 

x = R 3 x = ?f- = \ ( v^O^) + tan- 1 ^^4-) , (12) 
we can write 

dy 

Cy = -Ay + L 1 (x)x— + L (x)y, (13) 

where 

A _ d 2 N d 
dx 2 2 dx 

and L\(x) and Lo(^) are analytic on < x < 7r 2 i? 3 /16. While r runs over 
the interval [1, R], x runs over [0, xr], where xr := R 3 C, R /A(< tt 2 R 3 /16). The 
Dirichlet condition at the regular boundary is y\ x=XR = 0. 
Remark: Since x = i? 3 £ 2 /4, we have 

d 2 N d 1 / d 2 N - 1 d \ 
~ X dx^ + ~2dx~ ~R3\de + ~l~dj)' 
Thus A is the radial part of the Laplacian in the iV-dimensional Euclidean 
space M. N provided that N is an integer. But we do not assume that N is an 
integer in this study. 
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2.4 

Let us fix an eigenvalue A = X n and an associated eigenfunction <3>(x) of C 
Then 

yi(t,x) = sin y/\{t-t )${x) (14) 
is a time-periodic solution of the linearized problem 

d 2 y 



dt< 



+ £y = 0, y\ x =x R = 



Moreover we claim that <&(x) is an analytic function of \x\ < tt 2 R 3 /16. To 
verify it, we use the following 

Lemma 1 We consider the differential equation 

X dx>- + h{x) Tx +c{x)y = ^ 

where 

b(x) = (3 + [x]i, c(x) = [x] , 
and we assume that (3 > 2. Then 1) there is a solution y± of the form 

2/i = l+ [x] u 

and 2) there is a solution y 2 such that 

y 2 = x-^\l + [x] 1 ) 

provided that ft N or 

y 2 = x~ p+1 {\ + [x]x) + hy x \ogx 

provided that /3 EN. Here h is a constant which can vanish in some cases. 

For a proof, see PQ, Chapter 4. Applying this lemma with (3 = N/2 to 
the equation 



dx 2 



^-L 1 (x)x)^ + (X-L (x))y = 0, 



we get the assertion, since y 2 ~ x 2 cannot belong to X — L [x 2 dx) 
for N > 4, even if N — 4, which was the exceptional case in the preceding 
discussion of the limit point type. 
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3 Statement of the main result 



3.1 

We rewrite the equation (5) by using the linearized part C defined by (7) as 
d 2 y ( , „ ( dy\\ / dy 



dt 2 



+ {^ + Gi (y, r-£))Cy + Gn (r, y, r£) = 0, (15) 



where 



Gj(y,v) = {l + y) 2 (l + -d v G 2 (y,v))-l, 



7 

P ~ , s 11 



Gn(r, y, v) = —G II0 (y, v) + n G m (i/, v), 

pr z 7 — 1 r 6 

G no (y,v) = (l + y) 2 (3d v G 2 -d y G 2 )v, 

G in (y,v) = (1 + y) 2 (-(^G 2 )((4-3 7 )y- iv) + G 2 ) - H + Ay{\ + yf. 



Here 



G 2 := G- 1 (3y + v) = [y,v} 2 , 
d v G 2 := _G 2 = _- 7 =[ J/ , T ,] 1 

:= |-G 2 = ^-3 7 =| J/ , V ] 1 



We have fixed a solution of the linearized equation y tt + Cy = (see 
(14)), and we seek a solution y of (5) or (15) of the form 

y = ey 1 + ew, 

where e is a small positive parameter. Then the equation which governs w 
turns out to be 

d 2 w ( dw \ dw 

+ (1 + ea(t,r,w,r— ,e) jCw + eb(t,r,w,r— ,e) = ec(t,r,e), (16) 

where 

a(t,r,w,fi,e) = + u>)> e(«i + ft)), 

&(*, r, u>, ft, e) = e~ 1 G I {e{y 1 + w), e( Vl + ft))A/i + e- 2 G H (e( yi + w), e(ui + ft)) 
-e^G^eyi, ev 1 )Cyi - e~ 2 G n (eyi, e^i), 
c(f,r,e) = e -1 G/(eyi,eui)£j/i + e~ 2 Gjj(ej/i,eui). 
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Here V\ stands for rdyi/dr. 

The property of the trouble which comes from the free boundary touching 
the vacuum is same as that of the case of one-dimensional motions studied in 
[3] . Namely, because of loss of regularities at r = R, the usual iteration does 
not work for the equation (16), which is not quasi-linear despite appearances. 
So, we should employ the Nash-Moser theorem. 

Then the main result of this study can be stated as follows: 

Theorem 1 For any T > 0, there is a sufficiently small positive Eq(T) 
such that, for < e < £q(T), there is a solution w of (16) such that 



or a solution y G C°°([0,T] x [1, R]) of (5) or (15) of the form 

y(t,r) = e yi (t,r) + 0(e 2 ), 

or a motion which can be expressed by the Lagrangian coordinate as 

r(t, m) = f (m)(l + eyi(t, f(m)) + 0(e 2 )) 

forO <t < T, < m < M. 

Remark. The corresponding density distribution p = p(t,r), where r is 
the original Euler coordinate, satisfies 

p{t, r) > for 1 < r < R F (t), p(t, r) = for R F (t) < r, 



w G C°°([0,T] x [1, R}) and 




where 



R F (t) :=r(t,M) = J R + e J Rsinv / A(t-t )$(0) + 0(e 2 ). 
Since y(t, r) is smooth on 1 < r < R, we have 



pit, r) = C(t)(R F (t) - r)^i(l + 0{R F (t) - r)) 



as r — > R F (t) — 0. Here C(t) is positive and smooth in t. 
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3.2 

Our task is to find the inverse image < P _1 (ec) of the nonlinear mapping ^3 
defined by 

= "a^ + (! + + £& " ( 17 ) 

Let us note that ^3(0) = 0. This task, which will be done by applying the 
Nash- Moser theorem, will require a certain property of the derivative of 

d 2 h dh 
D<#(w)h = — + (1 + ea x )Ch + ea 21 r— + ea 20 h, (18) 
ot z or 

where 

, dw . 

ai = a{t,r,w,r—,e), 

da db 
aio = ^— ^ w + ^— , 

° 21 = + 

duu 

Here f2 is the dummy of r— — . The following observation will play a crucial 

Or 

role in energy eatimates later. 

Lemma 2 We /iave 

7-P, x 9^+9/ x . 9// 4 7^ 2e(7- 1) /0T\ 2 \ 

where 

dy 

y = e{y l + w) ) v = r—, Y = y 1 + w. 

Proof: It is easy to se 

oai = (duGrfCY + e- l d v G n 

= {d v Gj) f-^(3Y + V)') + e-iJLdvGm + — 

\ pr J pr 2 7 — 1 r 6 

where 

[[/] = 7 (^G / )(3F + V) + 9 w Gj(-4y) + e-^Gm- 
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Using 



we can show that [U] 
sion. ■ 



d v G I = (l+y) 2 -dlG 2 , 

7 

0. Then a direct calculation leads us to the conclu- 



4 Proof of the main result 

The proof of the main result is almost same as the discussion developed in 
[3], Section 4. We use the variable x defined by (8), (12) from r, and this x 
plays the role of the variable z in [3]. We note that 



x 
d_ 

Or 



R 2 (R-r) + lR-r] 2 , 



-R 2 {i + NO) 



d_ 

dx 



Therefore a function of 1 < r < R which is infinetely many times contin- 
uously differentiable is also so as a function of < x < xr. Therefore we 
apply the Nash-Moser theorem ([2], p. 171, III. 1.1.1) to the space of functions 
of t and x: 



<£ := C°°([0,T] x [0,x R ]), 
(S := {w G (£ | w\ t=0 = w t \ t =o = 0, w\ x=Xl 



0}. 



The nonlinear mapping ^5 maps <£q into (£. 

On <£(d <£q) we consider the equivalent gradings of norms: 



\y 



|(oo) 



sup \\{-^) 3{ - A)ky 

0<j+k<n 11 v Cft / 



Ml<? 



0<j+k<n 



T 



df- 



L^([0,T]x[0,x R ]) 
2 n 1/2 

it) 

x 



) J (-A) k y( x dt) 



Here y G (£, n G N, and X = L 2 ([0, x_r], x ^ Mx). Then (£ and it closed 
subspace (Eo are tame spaces. The mapping *p is tame thanks to the following 



Proposition 3 For any meNw have the formula 



A m Dw(x) 



JV 

X 2 " 



A m+1 w(x')(x')% 



+m-l 
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and 

1 

Here D = d/dx 



&) m D k w\\ L ~ < k _ lN IK-A^IUo 

n i=0 (f 



Therefore the problem is concentrated to energy estimates in the space 
X of the solution and its higher derivatives of the linear equation 

D%(w)h = g, 

when w is fixed in (£ - Let us investigate the structure of the linear operator 
First we note that 

J- = --- = — 1+ 4 ■ 

p r R R 6 

Therefore it follows from Lemma 2 that there exists a smooth function a(t, x) 
such that 

^a 2 ir— = ea{t,x)x—. 

Let us put 

&i := (1 + £ai)Li(rr) + ea, 
6 := (1 + eai)Lo(x) + ea 2 o, 

taking into account the observation in Section 2.3, (13). 
Then we can write 

d 2 h dh 
D^(w)h = — - (1 + ecu) Ah + &i(*> a;)a;— + b (t, x)h. 

We note that b\, b depends only on w, dw/dx, d 2 w/dx 2 . Then we can claim 
that the solution of D^(w)h = g enjoys the energy estimate 

\\h t \\x + \\Dh\\ x + \\h\\ x < C sup \\g(t)\\ x , 

0<t<T 

where D = y/xd/dx and C depends only on 

A = ||e:ai jt ||x,oo + y/2\\eDai + 6i ||x,°° 
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and 

B = H&oIIl 00 - 

A proof can be found in [3]. Similarly the higher derivatives H m = A m h 
satisfy the equations 

(H m ) tt - (1 + ea x ) AH m + b^ ] DH m + b { ™ ] H m = F m , 



where D = xd/dx and 



m— 1 

F m = A m g + J2(c i y ) DH j + c^H j ). 

j=0 

The recurrence formula to determine the coefficients b^ , b^ is the same 
as in |3j. But the initial condition is, for the present case, 

= &!, ^ = bo. 

The solution is 
&5 m) = -2mD(ea 1 ) + b 1 , 

b Q m) = (-m A - m(m - + 2D)D)(ea 1 ) + m(l + 2D)b x + b 

= -m((2m - 1)A + (m - 1)(1 - N)D)(ea 1 ) + m(l + 2D)b x + 6 - 

The coefficients 7 q™' ) = c^l_ k , 7 ^ = c^_ fc are determined by starting from 

7 {i ) = (-A + iV-2)6 1 , 7 £ ) = -A6 d , 
through the recurrence formula 

7ir 1} = 7!? + (A - (N - 2)D)^L + for 2 < k 

= # + 4m 2 (A + ^l^DjDCeoi) + 

+ (-(4m + 1)A + {2mN - 6m + N - 2)D)b t - 2Db 

= 7^ + (1 + 2^)7^ + A 7 ^ for 2<fc 

7 Sr +l) = 7 Sr } - A^ m) + (1 + 2D) 7 £r ) = 

= 7 <™ } + mA((2m - 1)A + (m - 1)(1 - N)D)(ea 1 ) + 
-m(3 + 2 J D)A6 1 - A6 + (1 + 2L>) 7 {™ ) 
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Thus if we put 

3 - N 

A = (Dai, Aai, (A + — - — D)Dai, Dh, Ab h Db , Ab ), 
we can show that 

ll7ir ) Ho+ \\l^\\o<C\\A\\ k . 

Here 

\\y\\e= sup sup \\(-Ayy\\ L oo. 
o<t<To<j<e 

This implies the estimate 

\\DH m \\ x + \\H m \\ x < C(\\g\\ m + \\A\\ m ), 
where C depends only on 

A m = \\ea ht \\ L oo + \Z2\\eDax + &i |U°o, 

which do not contain the higher order derivatives of w. This is the required 
tame estimate of the inverse of D^(w). 

The proof performed in [3] can be applied directly to the present case. 
This completes the proof of main result. 
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